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Abstract 

We present 10' x 50' scan maps around an M supergiant a Ori in at 65, 90, 140 and 160/im obtained 
with the AKARI Infrared Astronomy Satellite. Higher spatial resolution data with the exact analytic 
solution permit us to fit the de-projected shape of the stellar wind bow shock around a Ori to have the 
stand-off distance of 4.'8, position angle of 55° and inclination angle of 56°. The shape of the bow shock 

— 1/2 

suggests that the peculiar velocity of a Ori with respect to the local medium is v* — 40 n H ' , where hh is 
the hydrogen nucleus density at a Ori. We find that the local medium is of tih = 1-5 to 1.9 cm~ 3 and the 
velocity of the local flow is at 11 km s _1 by using the most recent astrometric solutions for a Ori under the 
assumption that the local medium is moving away from the Orion OB 1 association. AKARI images may 
also reveal a vortex ring due to instabilities on the surface of the bow shock as demonstrated by numerical 
models. This research exemplifies the potential of AKARI All-Sky data as well as follow-up observations 
with Herschel Space Telescope and Stratospheric Observatory for Infrared Astronomy for this avenue of 
research in revealing the nature of interaction between the stellar wind and interstellar medium. 

Key words: stars: individual (a Ori) — stars: kinematics — stars: mass-loss — stars: supergiants — 
ISM: structure 



1. Introduction 

Alpha Orionis (Betelgeuse, HD 39801; hereafter a Ori) 
is a supergiant of spectral type M2 lab: and is a very bright 
far-infrared (far-IR) source. Some extended structure in 
the far-IR was recognized around a Ori first by Stencel, 
Pesce, & Hagen Bauer (1988) in the All-Sky Survey data 
of the Infrared Astronomical Satellite (IRAS). Using an 
improved image reconstruction technique, Noriega- Crespo 
et al. (1997) produced images of better quality for a Ori 
and identified the extended structure around the star as an 
arc-shaped structure at the interface between the stellar 
wind and the local interstellar medium (ISM). 

Such structures can be formed around mass-losing stars 
at which the ram pressure of the ambient ISM balances 
with that of the stellar wind. These stellar wind bow 
shock arcs had already been found around hotter, more 
luminous OB stars in the IRAS All-Sky maps (Van Buren 
& McCray 1988). Most recently, Ueta et al. (2006) dis- 
covered a stellar wind bow shock arc around R Hya, an 

* Based on observations with AKARI, a JAXA project with the 
participation of ESA. 



asymptotic giant branch star (an evolved star of lower 
mass and lower rate of mass loss), using the Spitzer Space 
Telescope. 

Far-IR emission from these arcs is probably due mainly 
to thermal emission of cold dust components of the stellar 
wind bow shock whose temperature peaks at far-IR. There 
may be contribution from low-excitation atomic lines such 
as [Ol] 63 fim, [Ol] 145 /jm and [Cll] 158 /xm at these 
wavelengths. However, the exact emission mechanism of 
these far-IR bow shocks remains unclear, while an attempt 
to identify their spectroscopic nature is currently on-going 
with the Spitzer Space Telescope. Meanwhile, theoret- 
ical studies have been done both analytically (Wilkin 
1996; Wilkin 2000) and numerically (e.g. Mac Low et al. 
1991; Dgani, Van Buren, & Noriega-Crespo 1996; Blondin 
& Koerwer 1998; Wareing, Zijlstra, & O'Brien 2007) and 
the structure of the stellar wind bow shocks has been un- 
derstood reasonably well. 

The AKARI Infrared Astronomy Satellite (AKARI; 
Murakami et al. 2007) makes the All-Sky Survey in the 
far-IR for the first time in 25 years since IRAS at much 
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finer spatial resolution. 1 The original IRAS images of a 
Ori allowed only rough identification of the shape of the 
extended far-IR emission (Stcnccl, Pescc, & Hagen Bauer 
1988). The enhanced IRAS images of a Ori yielded just 
the mean radius and width of the arc with a rough esti- 
mate for the position angle of the apex (Noriega- Crespo 
et al. 1997). In this paper, we characterize the structure 
of the far-IR stellar wind bow shock around a Ori, using 
higher spatial-resolution AKARI images, and investigate 
the kinematics of the star, bow shock and ISM in the 
vicinity of the star by adopting the analytic solution of 
the bow shock (Wilkin 1996) and the most recent astro- 
metric solutions for the star (Harper, Brown, & Guinan 
2008). 

2. Observations and Data Reduction 

We observed a Ori in the four bands at 65, 90, 140 
and 160/xm using the Far-IR Surveyor (FIS; Kawada et 
al. 2007) on-board the AKARI satellite (Murakami et 
al. 2007) on 2006 September 21 as part of the MLHES 
("Excavating Mass Loss History in Extended Dust Shells 
of Evolved Stars") Mission Program (PI: I. Yamamura). 
Observations were made with the FIS01 (compact source 
photometry) scan mode, in which two strips of forward 
and backward scans were done with a 70" spacing at the 
15" s^ 1 scan rate, resulting in the sky coverage of 10' x 50' 
centered at the target. 

The FIS Slow-Scan Toolkit (Vcrdugo, Yamamura, & 
Pearson 2007 2 ; ver. 20070914) was used to reduce the 
data. We found that the quality of the resulting map was 
improved when we used a combination of the temporal 
median filter with the width of 200 s (or longer), tempo- 
ral boxcar filter with the width of 90 s, and sigma clipping 
threshold of 1.5. For the reduction of data in the short 
wavelength bands (SW bands; 65 and 90/im), the results 
were improved further when we performed flat-fielding us- 
ing the local "blank" sky data. Furthermore, we applied 
a custom reduction process to remove pixel-dependent re- 
sponse variation by subtracting the baseline "sky value" 
from each pixel determined by a linear least-squares fit to 
the off-source pixel values in the time-series data. 

The resulting maps are in 15 and 30" pixel -1 (default 
pixel scale) for the SW and LW (long wavelength; 140 and 
160/im) bands, respectively. The resulting 1 a sensitivities 
are 4.0, 4.1, 6.5 and 19.7 MJy sr" 1 while achieving, on 
average, five, eight, 15, and 10 sky coverages per pixel were 
at 65, 90, 140 and 160/mi, respectively. The sky emission 
(the component removed during the reduction) is found to 
be 32.3±4.6, 31.7±3.5, 56.1±5.8 and 81.4±9.8 MJy sr" 1 
at 65, 90, 140 and 160/im, respectively. Photometry was 
done following the latest calibration method to address 
the effects due to slow transient response of the Ge:Ga 
detectors (Shirahata et al. 2008). Image characteristics 
are summarized in Table 1. 

The measured sky emission values are consistent 

1 The effective beam size of AKARI (0.'5 to Of 9) is much smaller 
than that of IRAS (2.'0 to 5'.0). 

2 Available at http://www.ir.isas.jaxa.jp/AKARI/Obscrvation/ 



with the estimates obtained with the Spitzer Planning 
Observations Tool (SPOT) based on the COBE/DIRBE 
data (Reach 2000) in the wide bands at 90 and 140/im 
(28.7 and 52.1 MJy sr _1 , respectively), while about 50% 
larger in the narrow bands at 65 and 160/im (21.7 and 
53.8 MJy sr _1 , respectively). By design, the spatial scale 
of the measured sky emission corresponds at most to the 
scan length (50'), which is comparable to the spatial res- 
olution of the far-IR background data used in the SPOT 
background estimates (40' to 70'; Reach 2000). Hence, 
good agreement between the measured and estimated sky 
emission values in the wide bands is reasonable. The 50% 
discrepancy in the narrow bands, on the other hand, may 
indicate highly structured and variable [O l] line emission 
at 63 and 146/im in the background in the vicinity of a 
Ori, an oxygen-rich supergiant suffering from heavy mass 
loss. With our data, however, we are unable to either 
prove or disprove this possibility, unfortunately. Follow- 
up far-IR spectroscopic observations are indeed necessary. 

3. Results 

Figure 1 shows the background-subtracted false color 
AKARI/FIS scan maps of a Ori at 65/iin (N60; far left), 
90/im (WIDE-S; second from left), 140/mi (WIDE-L; sec- 
ond from right) and 160/im (N160; far right). Surface 
brightness in MJy sr -1 is indicated by the color scale that 
is specified in the wedge above each scan map. Due to the 
location of the detectors on the focal plane, the extent in 
the in-scan direction is different in each band, a Ori itself 
is clearly detected in all four bands. However, the 65 and 
90/mi band maps are affected by pixel cross-talk induced 
by the bright star (linear extension of the star oriented at 
64° east of north) while the 160/im band map is impacted 
by a ghost seen about 4' southeast of the star (Kawada et 
al. 2007). The extended central emission core is the topic 
of our forthcoming paper and will not be discussed here. 

Measured fluxes of the star (aperture defined by where 
surface brightness drops to the sky level) are 349 ±11, 
151 ±4, 32 ±1 and 49 ±2 Jy at 65, 90, 140 and 160/im, re- 
spectively. The flux values at 65 and 90/im are consistent 
with the previous IRAS measurements (Noriega-Crespo et 
al. 1997). These values are obtained by following the stan- 
dard method for point-source photometry including aper- 
ture correction elucidated in the Users Manual (Verdugo, 
Yamamura, & Pearson 2007), except for the last step in 
which the effects due to slow transient response of the 
Gc:Ga detectors is now addressed (Shirahata et al. 2008). 
The AKARI Ge:Ga detectors are known to underestimate 
the flux by roughly 60% due to slow transient response, 
and the latest calibration allows one to derive a correction 
factor via a power-law function of the total flux (source 
plus background) based on the calibration observations 
(Shirahata et al. 2008). Previously, the correction factors 
were assumed to be constants. 

Besides the central star, the arc and bar arc also clearly 
detected in a much finer spatial scale than in the previous 
IRAS maps (Noriega-Crespo et al. 1997), even though the 
scan width was not wide enough to capture these struc- 
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Table 1. Characteristics of the AKARI/FIS Map of a Ori 



Band 


A 

H 


AA 

[fim] 


Pix Scale 
[arcsec] 


Coverage 
[pixel -1 ] 


Csky 

[MJy sr" 1 ] 


[MJy sr" 1 ] 


oSPOT* 

lAsky 

[MJy sr -1 ] 


N60 


65 


22 


15 


5 


4.0 


32.3 ±4.6 


21.7 


WIDE-S 


90 


38 


15 


8 


4.1 


31.7±3.5 


28.7 


WIDE-L 


140 


■52 


30 


15 


6.5 


56.1 ±5.8 


52.1 


N160 


160 


34 


30 


10 


19.7 


81.4±9.8 


53.4 



* Estimate by SPOT (Reach 2000). 



tures in their entirety. Additional scans to cover the whole 
arc were scheduled on 2007 September 21. Unfortunately, 
cryogenic liquid Helium boiled off on 2007 August 26 and 
we were unable to obtain the additional data. While de- 
tailed analysis of the structure of the arc has to wait until 
the All-Sky data become available, we can still make use 
of the data at hand. 

The arc is distinctly visible to the north and south- 
east of the star in the 65 and 90^m band maps and is 
marginally discernible in the 140/im band map. In the 
160/im band map, however, the arc is blended in with 
the background cirrus and is not cleanly distinguishable. 
The bar, on the other hand, can be seen in all four maps, 
while the structure becomes progressively less well-defined 
at longer wavelengths. In the SW band images, surface 
brightness of the brightest parts of the arc (northeast of 
the star; 20 to 25 MJy sr -1 ) is higher than that of the 
bar (~ 10 MJy sr -1 ). In the LW band images, however, 
surface brightness of the bar (15 to 20 MJy sr -1 ) is higher 
than that of the arc (5 to 10 MJy sr 1 ). 

Assuming that emission from the arc and bar detected 
in these maps is mainly due to thermal emission of 
optically-thin concentration of dust grains, we can esti- 
mate the dust temperature by fitting the observed sur- 
face brightnesses at these bands to a Planck curve (i.e., 
S v oc 7-0(1/0 / 'v)P 'B u (T dust ) , where t„ = t q (is /v)P is the op- 
tical depth at v power-law scaled from To at vq with an 
index f3, B v is the Planck function at v and Tdust is the 
dust temperature). 

Instead of treating /3 as a free parameter, we varied it 
between 1 and 2 to see how Xd us t and To would behave. We 
find 42 ±9 K near the brightest parts to 11 ± 1 K near the 
edge of the arc and 22 ± 5 K in the bar. For the optical 
depth, we find values on the order of 10~ 5 to 5 x 10 -4 
typically with a factor of ~ 2 uncertainty. The optical 
depth in the arc is generally about an order of magnitude 
lower than in the bar. Table 2 summarizes the observed 
characteristics of the arc, bar and background sky. 

4. The Stellar Wind Bow Shock 

4-1. Orientation of the Bow Shock Cone 

The far-IR arc around a Ori has been interpreted as 
the interface between the interstellar medium (ISM) and 
the circumstcllar envelope developed by the stellar wind 
(Stencel, Pesce, & Hagcn Bauer 1988; Noricga-Crespo et 
al. 1997) at which the ram pressure of the ambient ISM 
balances with that of the stellar wind, as has been found 



Table 2. Characteristics of the Arc and Bar Near a Ori 





A 




?dust 






Object 


[/im] 


[MJy sr" 1 ] 


[K] 


TO 




Arc 


65 


5 to 25 


11 to 42 


4 x 10 


-5 




90 


5 to 20 




2 x 10 


-5 




140 


5 to 10 




8 x 10 


-6 




160 


10 to 15 




7x 10 


-6 


Bar 


65 


~ 10 


17 to 27 


5 x 10 


-4 




90 


- 10 




3 x 10 


-4 




140 


~ 15 




1 x 10 


-4 




160 


-20 




1 x 10 


-4 


Sky 


65 


36.2 


25 to 32 


5 x 10 


-4 




90 


37.9 




3 x 10 


-4 




140 


59.7 




1 x 10 


-4 




160 


85.9 




1 x 10 


-4 



around OB stars (Van Buren & McCray 1988). The shape 
of such bow-shock arcs has been shown to follows the curve 
z = 7 ,2 /3i?o (where Rq is the stand-off distance between 
the star and the apex of the bow and r is the distance of 
the bow from the symmetric z axis of the bow) both nu- 
merically (Mac Low et al. 1991) and analytically (Wilkin 
1996). 

The shape of the arc around a Ori, however, appears to 
be much more circular. If one traces the brightness peak 
of the arc to fit an ellipse, such an ellipse turns out to have 
the semi-major axis length of 9.6' and the eccentricity of 
0.02 centered at 4.0' off the star at the position angle of 
56° east of north. Mac Low et al. (1991) have already 
shown that the apparent shape of the bow becomes more 
circular upon the consideration of the inclination angle. 
Since we have an exact analytic solution for stellar wind 
bow shocks developed by Wilkin (1996), namely, 

R($) = R esc 0y/3(l -6 cote) (1) 

where 9 is the polar angle from the apex of the bow, and 
the observed shape of the arc is simply a conic section of 
an axisymmetric bow shock cone, we can fit the brightness 
peak of the arc to this solution. 

For a given set of the stand-off distance and inclination 
and position angles of the bow (these angles define the 
3-D orientation of the bow), one can predict how the bow 
appears in the plane of the sky assuming that (i) far-IR 
emission is mainly due to thermal dust emission, (ii) the 
shell is optically thin to far-IR light and (iii) the column 
density of dust is the highest where the bow intersects 
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with the plane of the sky. Dust radiative transfer calcula- 
tions in the circumstellar shells have shown that the above 
assumption (iii) is generally valid as long as the assump- 
tion (ii) is valid (e.g. Ueta et al. 2001; Ueta & Meixner 
2003). 

Then, one can quantify the difference of the prediction 
from the data by computing the inverse of the sum of the 
squares of the differences between the distance from the 
star to the arc peak and that from the star to the predicted 
positions of the arc. Since this quantity represents the 
"correlation" between the data and prediction, it tends 
to be small if the prediction differs from what the data 
suggest. Thus, one can find the best-fit parameter set by 
locating the point in the parameter space at which this 
quantity becomes the largest. 

Through this method, we find the stand-off distance 
(de-projected) to be 4.'8±0.'1, the position angle to be 
55° ±2° and the inclination angle to be 56° ±4°. The 
uncertainties stem from those in determining the bright- 
ness peak of the arc by fitting the Gaussian to the surface 
brightness profile of the 90pm map, in which the arc is 
the most well-defined. Figure 2 shows the position of the 
apex of the best-fit bow shock cone projected in the plane 
of the sky. 

The position angle is consistent with what has been 
concluded from lower resolution IRAS maps (60° ± 10°; 
Noricga-Crespo et al. 1997) and is now more accu- 
rately determined with higher resolution AKARI maps. 
Adopting the heliocentric radial velocity of u ra d = +20.7 ± 
0.4 km s _1 (the mean of the measurements made by Jones 
(1928) and Sanford (1933) as used by Harper, Brown, & 
Guinan (2008)), we see that the bow shock cone is oriented 
into the plane of the sky, unlike the previous studies as- 
sumed to be close to edge-on (e.g. Noriega- Crespo et al. 
1997; Harper, Brown, & Guinan 2008). This is the first 
observational determination of the inclination angle of a 
stellar wind bow shock, owing to higher resolution AKARI 
maps. 

4-2. Space Velocity of a Ori 

The stand-off distance is determined by the ram pres- 
sure balance of the wind and ambient ISM. Therefore, 
starting from the pressure balance equation p w v^ = pv%, 
one can determine the stand-off distance for a star that 
loses mass via an isotropic stellar wind of velocity v w and 
mass-loss rate M while traveling through the ISM of a 
uniform density p at velocity of v» , the stand-off distance 
is 



Ro 



1 Mv w 
4npv% 



(2) 



(Wilkin 1996). For a Ori, we have fairly well-established 
estimates for M = 3.1 ± 1.3 x 10~ 6 M Q yr" 1 (Harper, 
Brown, & Lim 2001) and v w = 17± 1 km s _1 (Bernat et al. 
1979). Also, Harper, Brown, & Guinan (2008) have de- 
rived the new distance of 197 ±45 pc to a Ori based on as- 
tromctric solutions obtained by combining Hipparcos data 
with multi-epoch, multi-wavelengths VLA radio positions. 



Using these values, the de-projected stand-off distance Rq 
is 8.5 ±1.9 x 10 17 cm. By keeping the interstellar hydro- 
gen nucleus density at a Ori, nn, as a free parameter, the 
peculiar velocity of the star with respect to the ISM in 
the vicinity of a Ori is 



47r^Hm H ".Hi?o 



(40 ±9) n H 1/2 (km s" 



(3) 



where pn is the mean nucleus number per hydrogen nu- 
cleus for local medium (~ 1.4) and mn is the mass of 
hydrogen nucleus. 

Given the orientation of the bow shock cone, the pe- 
culiar velocity of the star with respect to the ISM in the 
vicinity a Ori can be decomposed into each of the equa- 
torial space- velocity components (the radial direction, the 
direction in right ascension corrected for declination and 
the direction in declination) as 



- a Ori 

Vp 

v a = 

S J a Ori ISM 

These values can be 



33 ±8 
18 ±5 
13 ±3 



-1/2 



(km s" 1 ) (4) 



converted to the Galactic space- 
velocity components ([/, V and W where they are pos- 
itive in the directions of the Galactic center, Galactic ro- 
tation and the North Galactic Pole, respectively) through 
a spherical trigonometric transformation (e.g. Johnson & 
Soderblom 1987) as 

a Ori 

-34 ±8 

1/2 



u 

V 

w 



a Ori ISM 



-34 ±8 
-10±4 
-17±5 



(km s 



(5) 



Since these values are based purely on the orientation of 
the bow shock cone, they represent the Galactic space- 
velocity components of a Ori with respect to the ISM in 
the vicinity of a Ori. Here, the superscript to the [C/,V,W] 
vector indicates what velocity it refers to while the sub- 
script refers to with respect to what the velocity is defined. 

4-3. ISM Flow in the Vicinity of a Ori 

The new astrometric solutions by Harper, Brown, 
& Guinan (2008) yield the heliocentric Galactic space- 
velocity components of a Ori as 

a Ori 



u 

V 

w 







-22 ±1 
-12 ±3 
+21 ±4 



(km s" 1 ) 



(6) 



Thus, by combining the above two sets of values we can 
compute the heliocentric Galactic space-velocity compo- 
nents of the ISM in the vicinity of a Ori as 

a Ori ISM 



u 

V 

w 



' -22 ±1 " 




' -34 ±8 " 


-12 ±3 




-10 ±4 


+21 ±4 




+17±5 



n H 1/2 (kms- 1 )^) 



Hence, the stellar wind bow shock around a Ori is a con- 
sequence of a mass-losing M supergiant moving in the ISM 
that flows in the direction specified by the above space- 
velocity vector. 
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This ISM flow around a Ori must originate from some- 
where in the vicinity of a Ori. The most probable 
source of this ISM flow in the vicinity of a Ori is un- 
doubtedly the Orion OBI association (Warren & Hesser 
1977). Assuming that all four sub-associations (OBla 
through OBld) contribute equally to generate the ISM 
flow in the vicinity of a Ori, one can define the flow vec- 
tor toward a Ori at [X,Y,Z] = [— 182, — 66, — 31pc] em- 
anating from the mean position of the associations at 
[X,Y,Z] = [-350,-170,-T33pc]. 

Then, one can search for values of tih that would align 
the ISM flow vector [U, V, W]g ° ri ISM with this flow from 
the Orion OBI association within the quoted uncertainties 
in Eq. 7. Possible values of ?i H turns out to be 1.5 to 
1.9 cm -3 . These values are higher than nn — 0.3cm -3 
estimated for material in front of the Orion OB association 
(Frisch, Scmbach, & York 1990). In the direction of a 
Ori, however, the Hi column density is estimated to be 
0.27 x 10 22 cm -2 (Kalberla et al. 2005), implying < Uh > = 
4.4 cm -3 for the distance of 197 pc. Thus, our estimate 
of tt-h may not be too large. 

The heliocentric Galactic space-velocity components of 
the ISM How in the vicinity of a Ori can be converted 
back to the equatorial space- velocity components through 
an inverse spherical trigonometric transformation as 



a Ori ISM 







21 -33 n 
23 - 18 n 
9-13 n 



-1/2 
H 

-1/2 
H 



-1/2 
11 



(km s" 1 ) (8) 



and for njj = 1.5 to 1.9 cm 

a Ori ISM 



-6 to - 3 
9 to 10 
-1 to 



(km s 



star to the wind velocity (v*/v w ) is 1.7 to 1.9. Thus, 
the bow shock may be prone to instability (Dgani, Van 
Buren, & Noriega-Crcspo 1996; Blondin & Koerwer 1998) 
and may even develop vortices (e.g. Wareing, Zijlstra, & 
O'Brien 2007). 

Figure 2 is a close-up of the northern part of the arc 
overlaid with a line that delineates the Wilkin (1996) an- 
alytic solution at the 55° position angle at the 47° incli- 
nation angle. In general the Wilkin curve predicts where 
the surface brightness peaks along the arc extremely well. 
However, as one follows the arc structure from the apex 
to the downstream direction, there is a discontinuity of 
surface brightness along the arc at the position angle ~ 0° 
and a local enhancement of surface brightness at around 
(0', 7') from the star. The position of this local brightness 
enhancement is somewhat interior to the Wilkin curve un- 
like other parts of the arc. 

This drop of surface brightness accompanied by a lo- 
cal enhancement off the Wilkin curve may be due to vor- 
tex shedding caused by instabilities in the bow surface 
(Blondin & Koerwer 1998; Wareing, Zijlstra, & O'Brien 
2007). In an isothermal bow shock, Blondin & Koerwer 
(1998) found that instabilities would manifest themselves 
as wiggles in the bow on the length scale of the stand-off 
distance. Figure 2 shows that reduction and enhancement 
of surface brightness occurs about the stand-off distance 
away from the apex. Qualitatively, the observed struc- 
ture of the bow shock resembles to that of the numerical 
models showing vortex shedding and the development of 
vortex rings. Therefore, more detailed numerical investi- 
gations into the development of instabilities on the surface 
of the bow shock appears worthwhile. 

(9) ^.5. The Bar ahead of the Bow Shock Cone 



Thus, the ISM around a Ori flows at about 11 km s _1 
into the position angle of ~ 95° out of the plane of the sky 
(toward us). 

Since the stellar wind bow shock is a consequence of 
the peculiar motion of the star and the ISM flow at the 
star, the apparent orientation of the bow shock would not 
necessarily yield information on both. When the direction 
of the proper motion is aligned with the direction of the 
apex of the bow, the ISM at the star is stationary with 
respect to the Sun or flows against the peculiar motion 
of the star. On the other hand, when the direction of 
the proper motion is not aligned with the direction of 
the apex of the bow, the ISM at the star probably flows 
obliquely with respect to the peculiar motion of the star. 
Therefore, direct comparison of the orientation of the bow 
with respect to the direction of the peculiar motion of the 
star would provide a reasonable diagnostic for the presence 
and direction of the ISM flow around the star. 

4-4- Substructure of the Bow Shock 

For the derived values of jih = 1.5 to 1.9 cm -3 , the 
corresponding peculiar velocity of the star with respect to 
the ISM in the vicinity of a Ori is 33 to 29 km s" 1 (Eq.3). 
This means that the ratio of the peculiar velocity of the 



The existence of a linear bar structure ahead of a Ori 
is intriguing. At the moment there is no evidence that 
indicates the bar being a (by)product of the interaction 
between the ISM and stellar wind from a Ori or otherwise: 
the origin of the bar remains unclear. Similarly, there is 
no evidence that suggests the bar is co-spatial with a Ori. 
However, our data (Table 2) show higher optical depths in 
the bar than in the arc and in the background. Moreover, 
the bar and background seem to have roughly the same 
optical depth values. 

Thus, the bar appears to be more like the background 
cirrus than the arc, and hence, the bar is probably not 
caused by the interaction between the stellar wind and 
ISM. Also, assuming that these structures are made up 
with similar matter and have similar density, the differ- 
ences in the optical depth suggests that the bar is ex- 
tended along the line of sight like a sheet. If the bar 
represents a local concentration of matter co-spaced with 
a Ori, the motion of the star and angular separation be- 
tween the star and bar imply a collision between the two 
in about 3000 yr. Even if this is the case, it is still not 
clear whether the bar is caused by this ISM flow in this 
region or is caused by other external means (such as A 
Ori). 
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5. Summary 

AKARI/FIS scan maps around a Ori (10' x 50') at 65, 
90, 140 and 160/mi are presented. These images show the 
extended emission core and most of the circumstellar arc 
plus the northern bar structure at much higher spatial 
scale than in the previously obtained IRAS maps. Spatial 
resolution of the data is good enough to define the struc- 
ture of the arc to be fit with the exact analytic solution 
for stellar wind bow shocks, while the scan did not cover 
the arc in its entirety and there are some anomalies by 
pixel cross-talk in the SW band and by ghosting in the 
160/^m band due to the bright central star. 

Rather circular appearance of the arc (eccentricity 0.02) 
suggests that the stellar wind bow shock cone of 4.'8 de- 
projected stand-off distance is inclined at 56° with respect 
to the plane of the sky and oriented at 55° position angle 
(east of north). Adopting the distance of 197 pc, rate of 
mass loss at 3.1 x 10 -6 M Q yr _1 , and wind velocity of 
17km s" 1 , the peculiar velocity of the star with respect to 
the ISM at a Ori is found to be v* = 40 n H 1//2 km s -1 . 

By comparing the Galactic space-velocity components 
of the peculiar velocity of the star (based on the orien- 
tation of the bow shock cone) and that of the apparent 
motion of the star (based on the astrometric solutions), 
we derived the space- velocity components of the ISM flow 
around the star. Assuming that the ISM flow in the vicin- 
ity of a Ori emanates from the Orion OB 1 association, 
we find the particle number density per hydrogen nucleus 
(rin) m the ISM around a Ori to be 1.5 to 1.9 cm -3 , which 
translates to a 11 km s _1 flow. 

Owing to higher spatial resolution of the data, we may 
be witnessing the development of a vortex ring along 
the surface of the bow due to instabilities. This re- 
search demonstrates that far-IR images of stellar wind 
bow shocks are excellent diagnostic tools to investigate 
the kinematics of the bow and the characteristics of the 
ISM in the vicinity of the star. These stellar wind bow 
shocks have been found not only around M supergiants 
but also around OB stars (Van Buren & McCray 1988) 
and an AGB star (Ueta et al. 2006). Therefore, this av- 
enue of research is going to flourish with the coming of the 
AKARI All-Sky Survey (Murakami et al. 2007) followed 
by new opportunities with Hcrschcl Space Telescope and 
Stratospheric Observatory for Infrared Astronomy, be- 
cause there will be a wealth of new far-IR data on the stel- 
lar wind bow shocks around variety of mass-losing stars. 

This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France. Ueta is grateful 
to Drs. R. E. Stcncel and G. M. Harper for illuminat- 
ing discussions. Ueta also acknowledges contribution by 
a summer student, A. Karska, in developing some of the 
custom data reduction algorithms. 
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Fig. 1. AKARI/FIS false-color maps of a Ori in the SW 
bands - N60 (Gbiira) and WIDE-S (90/im) at ^"pixel" 1 
scale - and in the LW bands - WIDE-L (140/xm) and N160 
(160/^m) at 30" pixel -1 scale - from left to right, respectively. 
Background emission has been subtracted by a combination of 
temporal filters during data reduction. RA and DEC offsets 
(with respect to the stellar peak) are given in arcminutes. The 
wedges at the top indicate the log scale of surface brightness 
in MJy sr — 1 . North is up, and east to the left. 
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Fig. 2. Close-up of the northern arc structure in the WIDE-S 
(90/im) map. RA and DEC offsets (with respect to the stellar 
peak) are given in arcminutes. The wedges at the top indicate 
the log scale of surface brightness in MJy sr — 1 . North is up, 
and east to the left. The white line is the Wilkin (1996) 
curve that delineates the analytic solution at the 55° position 
angle at the 47° inclination angle. The black star and triangle 
indicate the position of a Ori and of the apex of the best-fit 
bow shock cone projected in the plane of the sky, respectively. 



